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ABSTRACT 

We present a new derivation of the equations governing the oscillations of slowly ro- 
tating relativistic stars. Previous investigations have been mostly carried out in the 
Regge- Wheeler gauge. However, in this gauge the process of linearizing the Einstein 
field equations leads to perturbation equations which as such cannot be used to per- 
form numerical time evolutions. It is only through the tedious process of combining 
and rearranging the perturbation variables in a clever way that the system can be 
cast into a set of hyperbolic first order equations, which is then well suited for the 
numerical integration. The equations remain quite lengthy, and we therefore rederive 
the perturbation equations in a different gauge, which has been first proposed by 
Battiston et al. (1970). Using the ADM formalism, one is immediately lead to a first 
order hyperbolic evolution system, which is remarkably simple and can be numerically 
integrated without many further manipulations. Moreover, the symmetry between the 
polar and the axial equations becomes directly apparent. 

Key Vifords: relativity - methods: numerical - stars: neutron - stars: oscillations - 
stars: rotation 



1 INTRODUCTION 

The theory of non-radial perturbations of relativistic stars has been a field of intensive study for more than three decades, 
beginning with the pioneering paper of Thorne & Campolattaro in 1967. Their work was in turn based on previous studies of 
black hole perturbations initiated by Regge & Wheeler in 1957. Because of the quite involved and tedious computations of the 
perturbed field equations, the main focus remained on non-rotating neutron stars, although the foundations for computing 
rotating relativistic stellar models have already been laid by Hartle in 1967. Only in the early 90s, work began to shift to 
the perturbations of rotating relativistic neutron stars. In most works, the slow-rotation approximation was used to tackle 
the problem. Chandrasekhar & Ferrari (1991) studied the axisymmetric perturbations, where they established the coupling 
between the polar and axial modes induced by the rotation (polar or even parity modes are characterised by a sign change 
under parity transformation according to ( — 1)' while the axial ones change as (—1)'"''^). Shortly later, Kojima (1992) presented 
the first complete derivation of the coupled polar and axial perturbation equations. These equations were the starting point for 
investigations of relativistic rotational effects on stellar oscillations and associated instabilities. Most of the work considered 
the simpler task of solving the perturbation equations in the frequency domain and, as a result, calculations involved the 
determination of the eigenfrequencies rather than the solution of the time dependent equations. A somewhat different approach 
based on a Lagrangian description was used by Lockitch et al. (2001) with the focus on the computation of rotationally induced 
inertial modes. 

Following the more than 40 years old tradition, it was quite common to work in the Regge- Wheeler gauge, although some 
groups were using different gauges or the gauge invariant formulation of Moncrief (1974). In particular, a gauge introduced 
by Battiston, Cazzola & Lucaroni in 1971 in a series of papers to study the stability properties of non-radial oscillations in 
relativistic non-rotating stars, has not received proper attention (Battiston, Cazzola & Lucaroni 1971, Cazzola & Lucaroni 
1972, 1974, 1978; Cazzola, Lucaroni & Semezato 1978). 

Since the perturbation equations of non-rotating stars are fairly simple, there is no real advantage of one gauge over the 
other. Note, however, that there was a long standing puzzle why in the Regge- Wheeler gauge the equations could be reduced 
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to a fourth order system, whereas in the diagonal gauge used by Chandrasekhar & Ferrari (1991), which was previously used 
by Chandrasekhar (1983) in the context of black hole perturbations, only a fifth order system could be obtained. The latter 
system gave rise to an additional divergent solution, which had to be rejected on physical grounds. This discrepancy was 
finally solved by Price & Ipser (1991), who showed that the diagonal gauge was not completely fixed and still possessed one 
degree of gauge freedom, giving rise to the additional spurious solution. 

As computer power has been enormously increasing within the last decade, the problem of evolving the fully non-linear 
3D Einstein equations in the time domain has finally become into reach of feasibility. Various groups around the world are 
now building robust codes that perform the time evolution of both single neutron stars or binaries in full 3D (see for instance 
Font et al. 2001, Shibata & Uryu 2001, and the review of Stergioulas 1998). Nevertheless, there is still considerable work in 
progress within perturbative approaches, which can provide us with a deeper understanding of the physics of relativistic star 
perturbations. Moreover, any trustable non-linear code must be able to reproduce the results from perturbation theory. 

Time evolutions of the perturbative equations have been carried out first for the axial equations (Andersson & Kokkotas 
1996) and then for the polar equations using the Regge-Wheeler gauge by Allen et al. (1998) and Ruoff (2001a). Allen et 
al. (1998) managed to write down the evolution equations as two relatively simple wave equations for the metric perturbations 
and one wave equation for the fluid enthalpy perturbation inside the star. Ruoff (2001a) rederived these equations using the 
ADM formalism (Arnowitt et al. 1962). They were used to evolve and study initial data, representing the late stage of a 
binary neutron star head-on collision (Allen et al. 1999). 

Using the Arnowitt-Deser-Misner (ADM) formalism, the evolution equations for the axial perturbations of rotating stars 
have been brought into a suitable form for the numerical evolution by Ruoff & Kokkotas (2001a,b). Here, the resulting system 
of equations came out immediately as a first order system both in space and time, and it could be directly used for the 
numerical evolution without many further manipulations. In the non-rotating case, it is an easy task to transform the first 
order system into a single wave equation for just one metric variable. In the rotating case, however, this is not possible any 
more because of the rotational correction terms. 

When looking at the set of polar equations derived by Kojima (1992) it appears clearly that the presence of mixed spatial 
and time derivatives makes them not suitable for the numerical time integration. Nevertheless, using a number of successive 
manipulations and the introduction of many additional variables, we were able to cast the equations into a hyperbolic first 
order form. 

A more natural way to automatically obtain a set of equations, which is first order in time, is to use the ADM formalism. 
Ifowever, as we shall explain, even in that case the polar equations in the Regge-Wheeler gauge need to be further manipulated 
before they are suitable for a numerical integration. In general the ADM formalism yields a set of partial differential equations 
which are first order in time, but second order in space. For the numerical evolution, this is not ideal and one would rather 
prefer to have a pure first order system, or if possible a pure second order system, thus representing generalized wave equations. 
As we mentioned above, in the non-rotating case, it is easily possible to convert the perturbation equations into a set of wave 
equations. However, in the rotating case, this is not possible any more, even in the simple case when only axial perturbations 
are considered. To illustrate the problems associated with the Regge-Wheeler gauge, let us recall Einstein's (unperturbed) 
evolution equations written in the ADM formalism: 

(9t-£/5)7.. = -2aK.j , (1) 
{dt-C0)K,j = -a,,,+a[R,j+K''\K,j-2K,kK'',-A-K{2T,j-T\-i,j)\ , (2) 

with a denoting the lapse function, /3* the shift vector, Cp the Lie-derivative with respect to /3\ 7ij the metric of a space-like 
three dimensional hypersurface with Ricci tensor Rij, and Kij its extrinsic curvature. It is obvious that the only second order 
spatial derivatives are didjct and didjjki with the latter coming from the Ricci tensor Rij. This is still true for the linearized 
version of Eqs. (|^) and (^). 

In the Regge-Wheeler gauge, we have a non-vanishing perturbation of the lapse a and of the diagonal components of the 
spatial perturbations hij. Using the notation of Ruoff (2001a), the polar perturbations can be written as 

a ~ ^5r(t,r)y„„(e,0) , (3) 

Slr{t,r) \ 

n"'it,r) \Yi^ie,^). (4) 

sin2 6irj'"(t,r) / 

The perturbation equations coming from Eq. (^) contain second order r-derivatives of Si and T2. Note that they do not contain 
second order derivatives of S3, because only the angular components of the metric get differentiated twice with respect to 
r. In the axial case there is only one perturbation function for the angular metric components, but it is set to zero in the 
Regge-Wheeler gauge. Therefore the ADM formalism immediately yields a first order system. 

The polar equations, in contrast, can be cast into a fully first order system only through the introduction of some auxiliary 
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variables. In the non-rotating case, this is a fairly easy task, but for the rotating case, it turns out to become considerably 
more complicated. This is because the simple proportionality of Si and S3, which holds in the non-rotating case and which 
makes is easy to replace Si by S3, does not do so any more in the rotating case. Instead this relation involves rotational 
correction terms, and the replacement of 5*1 by 5*3 would lead to a considerable inflation of the equations. 

Consequently, instead of manipulating the perturbation equations in the Regge-Wheeler gauge, we look for a gauge in 
which the perturbation equations, by construction, do not show any second order spatial derivative. We have seen that the 
second derivatives come from the angular terms in the spatial metric and the lapse function. It seems therefore natural to 
set these components to zero. For the axial case this is already realized in the Regge-Wheeler gauge. It is only for the polar 
perturbations that we need a different gauge. From the seven polar components of the metric, Regge & Wheeler (1957) chose 
to set the components Vi, V3 and Ti, which represent in the notation of Ruoff (2001a) the polar angular vector perturbations 
and one of the angular tensor perturbations, to zero. We now choose a different set, namely we set the angular terms in the 
spatial metric Ti , T2 together with the lapse 5*1 to zero and retain Vi and V3 . With this choice we expect the ADM formalism 
to provide us with an evolution system without any second r-derivatives. 

We should mention again that this gauge has actually been introduced thirty years ago by Battiston, Cazzola & Lucaroni 
(1971) to derive the perturbation equations for non-radial oscillations of non-rotating neutron stars and to investigate in a 
subsequent series of papers their stability properties (Cazzola & Lucaroni 1972, 1974, 1978; Cazzola, Lucaroni & Semezato 
1978a, 1978b). Particularly relevant is the first paper of the series, in which they proved the uniqueness of this gauge, hereafter 
referred to as the BCL gauge, and where they also showed the relation with the Regge-Wheeler gauge. 

In Section 2, we will use the ADM formalism to derive the time dependent perturbation equations for slowly rotating 
relativistic stars in the BCL gauge. Section 3 contains a brief discussion of the non-rotating limit and conclusions will be 
given in Section 4. In the Appendix, we present the perturbation equations as they follow directly from Einstein's equations 
in a form similar to the equations in the Regge-Wheeler gauge given by Kojima (1992). Throughout the paper, we adopt the 
metric signature ( — h ++), and we use geometrical units with c — G = 1. Derivatives with respect to the radial coordinate r 
are denoted by a prime, while derivatives with respect to time t are denoted by an over-dot. Greek indices run from to 3, 
Latin indices from 1 to 3. 



2 THE PERTURBATION EQUATIONS IN THE ADM FORMALISM 

The metric describing a slowly rotating neutron star reads in spherical coordinates {t, r, 6, 










(5) 



where v, A and the "frame dragging" u; are functions of the radial coordinate r only. With the neutron star matter described 
by a perfect fluid with pressure p, energy density e, and 4-velocity 



(7" = (e-',0,0,ne' 



(6) 



the Einstein equations together with an equation of state p — p{e) yield the well known TOV equations plus an extra equation 
for the frame dragging. To linear order, it is given by 



7" — ^4nre'^^{p + e) ^ uj' — IGvre^^ (p + e) zu = , 



where 

■cu := Q — u) 



(7) 



(8) 



represents the angular velocity of the fluid relative to the local inertial frame. In the language of the ADM formalism, we have 
to express the above background metric (^ in terms of lapse , (covariant) shift and the 3-metric, which we denote by A, Bi 
and "jij, respectively. Explicitly, we have 



A = yjB'B, - goo = e" -f 0{lu^) 
Bi = (0, 0, -Ljr^ sin^ , 



7ij = 












(9) 
(10) 

(11) 
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The extrinsic curvature of the space-like hypersurface described by 7ij can be computed by 



yielding the only non-vanishing components 
1 



K 



K31 



(12) 



(13) 



The perturbations of the background lapse A, shift Bi, 3- metric 7^^, extrinsic curvature Kij, 4- velocity Ui, energy density e 
and pressure p will be denoted by a, Pi, hij, kij, Ui, 5e and Sp, respectively. The twelve evolution equations for hij and kij are 
obtained by linearizing the non-linear ADM equations (jl|) and (|^). Working in the slow- rotation approximation, we keep only 
terms up to order Q. (or uj). The background quantities B'' and Kij are first order in f2, hence we can neglect any products 
thereof. Furthermore, it is A' = 'y^-' Kij — 0. These circumstances lead to cancellations of various terms and the perturbations 
equations reduce to: 



dthij 
dtkij 



where 
k 

sr'^ij 

5R,i 



d,Pj + - 2 {Ak,, + K,ja + + Bfc(5r\,) , 

a [Rij + A-k{p — e)"fij] - didja + F^jdka + SF^jdkA 

+A [SR^j + K^jk - 2 (Ki'^kjk + K/'hk) + 47r ((p - e)hij + 7»j {5p - Se) - 2(p + e) {u^Suj + UjSur))] 



+B''dkkij + IdkKij - Ki'dj^ki 



K/^^-/kl) /?'■' + hkdjB'' + kjkd^B'' + K,''djPk + if/ftA- 



7'^ k- 



^3 I 
km 



7 [dihm] + djh 

n 



2r',;,/t, 



ij ^''Im j 1 



d,5T\k + t\ 



5v\k + r\k5r\ 



(14) 

(15) 

(16) 
(17) 
(18) 



To obtain a closed set of evolution equations, we will also use the four evolution equations the fluid perturbations following 
from the linearized conservation law iSr^% = 0. Last but not least we need the four linearized constraint equations, which 
serve to construct physically valid initial data and to monitor the accuracy of the numerical evolution; 



h''R, 



2K''k, 



167r {Se + 26"" (p + £)(«- 0^)5^3) 



(19) 



- Stt [{p + e)5ui + Ui {5p + Se)] = "f^^ {dikjk - djhk - T\kkji + V^khi ~ 5T\kKji + SV'jkKu) 

-h'" {d,K,k - djK.k ~ T\kKji + V'^kK^l) . (20) 

We assume the oscillations to be adiabatic, thus the relation between the Eulerian pressure perturbation 5p and density 
perturbation Se is given by 



Sp = 



Tip 
p + e 



Se + p'C[^-l 



where Fi represents the adiabatic index of the perturbed configuration, F is the background adiabatic index 
p + e dp 



P 



de ' 



(21) 



(22) 



and is the radial component of the fluid displacement vector The latter is related to the (covariant) 4- velocity pertur- 
bations Sun as follows 



C 1^ J, 1 A dS^" 1 K \i 

For the r component, this gives us 
{dt + 09^) C = e~^^ {e'Sur — Pr - ^hr,/,) . 



(23) 



(24) 



To proceed further, we expand the complete set of perturbations variables into spherical harmonics Yi^ = Yim{0, (j)). This will 
enable us to eliminate the angular dependence and obtain a set of equations for the coefficients, which now only depend on t 
and r. It is only then that we can finally choose our gauge. In principle, choosing the gauge amounts to providing prescriptions 
for lapse a and shift Pi. Those so-called slicing conditions determine how the space- like 3- metric foliates the 4-dimensional 
spacetime. In perturbation theory, the gauge can be used to set some of the ten metric perturbations to zero. We could, for 
instance, set a = /3i = 0, and we would be left with only the six components hij. Note, that setting a to zero is possible, since 
this is only the perturbation of the background lapse A, and the latter does not vanish. 

However, our actual goal is to set some of the spatial perturbation components hij to zero, in particular the angular 
components hab with a,b = {0,<f)}. In principle we can prescribe the values of hij only once for the initial data, and not 
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throughout the evolution. The only possible way to keep hat zero throughout the evolution is to choose our gauge such that 
the evolution equations for hab become trivial, i.e. we have to have 

dthab = 0, a,b^{e,4>}. (25) 
We will see that this requirement leads to a unique algebraic condition for the shift vector With the definitions 

Xi„, ~ 2 {de^ cot 6) d4,Yi^, (26) 

Wim := ( ale - cot edg - ]Yi^ ^ { 2die + /(/ + 1) ) , (27) 



^ Yira = {2dle + l{l + 1)) , 



we now expand the metric as follows. For the polar part we choose (symmetric components are denoted by an asterisk) 

a = , (28) 



<polar 



F polar 



E 



e bo , 



e A3 






* 








* 









and the axial part is 



E 

l.m 



V smt* 



E 

l.m 



-vr 



Vrsm6de\Yi 



Yw. 



(29) 
(30) 

(31) 
(32) 



For the extrinsic curvature we have all six components 



polar 



l.m 



1 

- -e X 

e Ki Yir. 

★ 



AKir) Ylrr, + KtWlr, 
-"■5 Aim 



e K2 04,Yim 

-K5 Aim 



1 _ 

-e 

2 



■E 

l,m 



■e Kg . ^ 



Im A;„ 



sin 6 



K'fT sine W,m 
KtsinOXim 



(33) 



(34) 



y * sine Wim 

Herein and throughout the whole paper, we use the shorthand notation 

A:=/(; + l). (35) 

We should note that the somewhat peculiar looking expansions for the coefficient Ki^"^ can actually be written as 

Wim - tCYim = 2dleYim , (36) 
-sm'eiyVim^KYim) = 2{coseiimede + dl^)Yim, (37) 

which are essentially the diagonal terms of the Regge- Wheeler tensor harmonic (c.f. Eq. (20) of Ruoff 2001a). However, 
we prefer to write them in terms of Wim. and Yim because it is only for these quantities that simple orthogonality relations 
apply. Furthermore, we have to mention that in the definition of the polar components of the extrinsic curvature, we differ 
from the notation of Ruoff (2001a), where the meaning of and is reversed (c.f. Eq. (24)). Also, the expansion for the 
axial perturbations is not exactly the same as in Ruoff & Kokkotas (2001a,b). 

In their original paper, Battiston et al. (1971) did not use the ADM formalism to fix the gauge, instead they defined their 
gauge by directly setting htt, hee, hg^ and h^^ to zero. Since the relation between htt and the lapse a is given by 

htt = 2Aa + 2B'/3i = 2e"a - 2uht^ , (38) 

it follows that in the rotating case, htt 7^ although the lapse a vanishes. In the non-rotating case l3i = and both a and 
htt vanish. If we insisted on keeping a vanishing htt also in the rotating case, we would obtain a non-vanishing lapse, giving 
us undesired second order spatial derivatives in the perturbation equations. 
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Finally the fluid perturbations are decomposed as 



Su; 



polar 



5e = 
5p — 

c = 



-e 2^[u\ ,1*2 oe,U2 a,^j Kim , 

-e > 0,-1/3 ^m^BdeXYir. 

P Ylm , 

l.m 



From Eq. (|l|), we have the relation 

Im _ (P + e) ^ rri 

For the sake of notational simplicity, we will from now on omit the indices I and m for the perturbation variables. With the 
above expansion, the evolution equations for hij read: 



(39) 
(40) 
(41) 
(42) 

(43) 
(44) 



{dt + imuj) SiYim 

J Ir 



dt{Vid4, + Vismede)Yi 



= (25'^ + 2A'S2-/Ci)y,„ 

+ 2uoe-'^^ {Vz ~ X'Vi) d^Yim + 2ue~'^^ {Vl - AV4) sinSaeY,™ , 

= (Vi~ -V, + e'^ (52 - K,)) deYr. ^ (vi ^ -V, - e^'K,) ^ 
\ r J \ r J smO 

— LJAV4 sin 6Yi,-a , 

= f Vl' - - Vl + e^^ (^ - ^4:2)) d4>Yim + (V2 - -V2 - e^'-Ks) smOdeYi 



(2S2 - ^Vi -K4 + ^A-s) Vim + 2Loe-^^ {V3d4.Yra + V4 sinedgY,^) 



= {Vi~ K5)Wi^ + {V2~ Kg) 



(45) 

(46) 
(47) 

(48) 

(49) 
(50) 



= ( Vl - J-^s) - {V2 - Ke) sin eWim ■ 

Still, in every equation a sum over all / and m is implied. From Eqs. (^ ) and (^^ we immediately obtain our condition for 
the shift components 



Vl = K,, 
V2 = Ke, 

and from Eq. (^) it follows after multiplication with F;^ and integration over the 2-sphere that 
S2 = ^K4 - uje-^^ {imVs + CfWi) , 



(51) 
(52) 

(53) 



where we have deflned the operator Cf^, which couples the equations of order I to the equations of order I + 1 and Z — 1, as 
■- ^^i'm' / Yil, sin edgYi,^,dn = (Z - l)(5i„yli_i,„ - (Z + 2)(3;+i„A;+i„ , (54) 



*-l ^Im 



with 

Qlm 



{I -m)(l + m) 



{21 - 1){21 + 1) 
Later, we will also need 



*-'2 ^Im 



and 



*-3 -^Im 



deY^mSvadYii^idn — —{l + l)Qi TTL-^l — lm + lQl + \mAl+\r, 



(55) 



(56) 



^^i'm' + Ytl, cos eYi,^,d^l + J Ytl, sin e^Yi,^,dQ 

I'm' \ S2 S2 

{I - 1)(/ + l)Oim^;-lm + l{l + 2)Qi + l„Ai + l„ . 



(57) 
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The operator C^^ can actually be expressed in terms of Cf^ and C^^: 



1 



(58) 



By making use of these relations we can eliminate the spherical harmonics and obtain the following simple set of evolution 
equations for the metric perturbations: 



{dt+imuj)S3 = Ki- Ki + X'K4-2uj'e'^^ {imV3 + jCfWi) , 

1 

— f 
2 



{dt + imuj) Va = Ki, - e^^Ki - -K5 + ^e^^Ki 



{dt + imLo) Vi = e^'^Ki - -Ke 



(59) 
(60) 
(61) 



In a similar way, we obtain the evolution equations for the six extrinsic curvature components, which are a little more lengthy: 



{dt + imuj) Ki 



{dt + imuj) K2 = 



_ y _ 

r r 



" - 1 , ^2A A 

J.2 '2r'^ 



S3 



+ STre'" {p + e) C^' [(cf ^ l) H + ^] - 2e 



2\ I 



im 



(62) 



+ -] Sz-^Vz 



.A- 1 



' [Ki-Ki + \'K4 - 4 — —K5 - 16nzu{p + e) {e^^ Ka + 2e^"ui) 



{dt + imuj) Kz 



{dt + imuj) K4, 



{dt + imuj) = e 



-^^i:fH{A-2)Ke) , 

^2.-2A A-2^^ ^ ^-2A^ /2.^^^ + (A - 2) £±^7^5 



— re •'-2 
2A ^ 



y - + X'Ki) - 167r7J7(p + e) (e'^^K^ + 2e^"ui) 



2A, 



+ STrre"" (p + e) CT' [(C? - l) H + 
+ r (£±^ - £±^) (cj'ATa + l&Tie^" vj{p + 6)^3) , 
(l/3'+(^'-A')V3-ie^^S3) 

+ |im Lj' (^^-^^4 — — IQ-Ke^" u^{p + e)u2 — £^^ (ijj'-R's + 167re^''ro(p + e)uz^ | , 



{dt + imw)Ka = e^-^^ (Vi + - A') 14) 



|im ([^'i^3 + IQ-Ke^" ^{p + e)u3) + J ^^-^4 — -^"2^ — l^-Ke'" ■uj{p + e)M2 J- 



(63) 



(64) 



(65) 



(66) 



(67) 



It is worthwhile to point out the symmetry between the polar and axial equations. Each pair V3 and V4, K2 and Kz, and K-o 
and K(, represent the polar and axial counterparts of a metric or extrinsic curvature perturbation. Thus, each associated pair 
of equations ( |60| ) and (1]), (H) and (|^), and (H) and (|^ has basically the same structure, with only the polar equations 
containing additional terms as there are more polar variables than axial ones. 

The last missing set of evolution equations is the one for the fluid quantities, coming from 5 (T'"'.,,) = and from Eq. (p^: 



{dt + vm^) H = 



+ rue 



I , ( n I \' , 2\ 2A A 

V r / 

,(V3'+0_A')V3 



K, A ,^ 1 ,^ 
+ -Ai - — ifs + --?f4 



H--Sz\\^C 



±^(i4'+(^-A')t4) 



e^" ^^iti + -A"4 + roe (imy3 + -Cf ^Vi) 



(ft + imS7) u\ = 



1 ip 



^^"ro (if^ - ^ifs) + (a;' + 2ro (v' - i 
"ro I^A's — - + ^cij' + 2td {v' ^ ) 



112 



(ft + imf7) M2 = // + ^ (im (27^2 - e-'" (A - 2) A5) + 2£±'u3 - e'^" L^^ ((A - 2) As)) - 



(68) 

(69) 
(70) 
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{dt + imfi) U3 = 2 — 



' ("3 + 



——Co A 
A 



{dt + imfi) — u' 
where 



e lii + -A4 + ii7e [imVa + V4) 



(71) 
(72) 



A 



vjC, < e 



A + (2Z.' - A' + ^ 



2A A 



+ e" 



2 r 



ti7 I 

r 



-2A . 1 -2l^ 

e Ml + -e A4 



(73) 



(74) 



In Eq. (pSj), we have defined the sound speed Cs as 
2 _ Ti dp 

The evolution equations comprise fourteen equations in total: four axial and ten polar ones. In the non-rotating case, they 
are equivalent to four wave equations, one for the axial and two for the polar metric perturbations plus one wave equation 
for the fluid variable H. The fluid equation for the axial velocity perturbation 113 vanishes in the non-rotating case, whereas 
equation (|7^) for the displacement variable ^ does so in the barotropic case. 

Finally we have the four constraint equations. The Hamiltonian constraint reads 

Syrr^e^^p = rS'3 - KVi + f 1 - 2rA' + )-e^^A) Sa + A f A' - i W3 



2 2A 

+ r e 



im (^^'^'e ^" K2 + IGirzu (p + e) U2^ + Cf^ {^'^'^ ^'^ K3 + Wtv-cj (p + e) u-j 



(75) 



and the three momentum constraints are 



K - -K's - A'l + ^'^K2 + 4(1+ ^'^O - i^'K, + '^rJS, 
r 2r r'' ' 4 



16nre " {p + e)u2 



IGnre^'^ {p + e)u3 = 



- (Syrr (p + e) + 2e"^^a;') {imVa + -Cf V4) , 
-rK'i + rKi + irv' - rX' - 2) K2 - -K5 + Ki + e' 



-2\ru) 
A 



(2imV3 - (A - 2) Ci^Vi) 



A 



1 -2A I r,l 

— e 
2 



u'S'a ~ 16^tu(p + e) {Sz + C:^ {{d + l)H~()) 



rK's + (ru - r\' - 2) K3 + ^^-^Ka + e"^^^ {2imV4 + (A - 2) CfWz) 



A^ 



\e-^^uo'S'z - l&nu,{p + e) (S3 + C'^ ((Cf + l) H ^ ()) 



A 

-2 1 1 r^2 



(76) 



(77) 



(78) 



The axial equations (|6l|), (|64[), (p7|), ( [ri[ ) and (^) without the coupling terms to the polar perturbations are equivalent 
to Eqs. (7)-(10) and (12) of Ruoff & Kokkotas (2001b). Note, however, that therein a slightly different definition of the 
perturbation variables has been chosen. 



3 THE NON-ROTATING LIMIT 



Although the non-rotating limit is well described by the wave equations given by Allen et al. (1998), it is instructive to 
consider it in the BCL gauge. This is obtained by setting f2 and ijj to zero in all the evolution equations (^9|)-([72[) and the 
constraints (^^~(^8|). As is well known, in this case the polar and axial parts of the equations completely decouple. For 
barotropic perturbations (Fi = F), the polar evolution equations can then be easily transformed into three wave equations 



for the rescaled metric variables 5* ; 



t^-^S-i and V 



e" Vjj/r and the rescaled fiuid variable H = e " H/r: 



9*2 



d S 2v-2\ 

+ e 



-I- 87re' 



+ 



A' 



9r* 



(^^'-A')+3 

(C?-l) (p'+(.'-i)p) + (cf)' 
A' 



,A 



A" 5 + 



4A 



V 



+ 2- 



V 



^ + 
r r 



(79) 



(80) 
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e ■< C, 



+ 



Qj.2 



Cl ( A'(:^ + A') +■ 

r 



dr 

2\ 



X" 



A 



A V V 



V + A') 



H 



+ e 



/A {c! ~l)V 



„2A 



- 1 



In Eqs. and (po|), r, is the well-known tortoise coordinate, which is related to r through 



d 



dr 



Furthermore, one can express the energy density p in terms of H, which in the barotropic case reduces to 



CI 



-H . 



(81) 



(82) 



(83) 



Although the equations in the first order form are quite simple, the above set of wave equations is more complicated than the 
equivalent set in the Regge-Wheeler gauge (Eqs. (14), (15) and (16) of Allen et al. 1998). This is particular so for the way in 
which the fluid variable p (or equivalently H) couples to the metric variable S, where the derivative of p enters. If the stellar 
model is based on a polytropic equation of state p = , then the behaviour of p at the stellar surface strongly depends 
on the polytropic index T. As discussed in Ruoff (2001a), p actually diverges for F > 2. In this case the metric quantity S 
would not even be C°, which can be troublesome for the numerical convergence. Although this could be a drawback of the 
BCL gauge, there is a clear advantage if one is interested in computing the gauge invariant Zerilli function Z in the exterior. 
Following Moncrief (1974), the definition of the Zerilli function is 



Z = 



with 



Afei 



'4e-4^fc2) 



r (A - 2) + 2M 



and 

I f BX — f r» 

fc2 = -e S . 



In terms of S and V this gives us 
■(S-AV^) 



(84) 



(85) 



(86) 



(87) 



r (A - 2) + 2M 

which is a simple algebraic relation in contrast to the relation in the Regge-Wheeler gauge, which includes a spatial derivative 
of one of the metric perturbations (see Eq. (20) of Allen et al. 1998, or Eq. (60) of Ruoff 2001a). In the Regge-Wheeler 
gauge, the two metric variables {S and F in the notation of Allen et al. 1998, and 5* and T in the notation of Ruoff 2001a) 
have different asymptotic behaviour at infinity, in particular one {F or T) is linearly growing with r. It is only through the 
delicate cancellation of the growing terms that the Zerilli function remains finite at infinity. However, this cancellation can only 
happen if both metric variables exactly satisfy the Hamiltonian constraint. Any (numerical) violation leads to an incomplete 
cancellation, and the Zerilli function starts to grow at large radii. This makes it very difficult in the numerical time evolution 
to extract the correct amount of gravitational radiation emitted from the neutron star. With the above relation (^) , we do 
not expect such difficulties to occur. 



4 CONCLUSIONS 

We have presented the derivation of the perturbation equations for slowly rotating relativistic stars using the BCL gauge, 
which has been first used by Battiston, Cazzola & Lucaroni in 1971. This gauge is defined by setting a, hee, hg^ and h^^ to 
zero. In the non-rotating case, the vanishing shift condition leads to a complete vanishing of htt, however, in the rotating case 
htt becomes non-zero (see also Appendix) . The advantage of the BCL gauge over the Regge-Wheeler gauge is that in the ADM 
formalism, the evolution equations a priori do not contain any second order spatial derivatives. Instead, one is immediately 
lead to a hyperbolic set of first order evolution equations, which can be directly used for the numerical time evolution without 
many further manipulations. Although it is in principle also possible to derive a hyperbolic set in the Regge-Wheeler gauge. 
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the procedure is rather tedious and requires the introduction of carefully chosen new variables in order to replace the second 
order or mixed derivatives. 

The perturbation equations for slowly rotating relativistic stars form a set of fourteen evolution equations plus four 
constraints. In the non-rotating barotropic case, it is possible to cast the equations into a system of four wave equations, three 
for the polar and axial metric perturbations and one for the polar fluid perturbation, as it is the case in the Regge-Wheeler 
gauge. Although these wave equations are not simpler than the corresponding ones in the Regge-Wheeler gauge, the first order 
system actually is. Maybe the main advantage is the simple algebraic relation of the metric variables to the Zerilli function. 
It was demonstrated by Ruoff (2001a) that the accurate numerical evaluation of the Zerilli function in the Regge-Wheeler 
gauge is somewhat difficult and requires high resolution because a small numerical violation of the Hamiltonian constraint 
can lead to very large errors in the Zerilli function. This should not be the case in the BCL gauge as the relation (|8^) does 
not involve any derivatives. 

A further advantage of these evolution equations is that the inclusion of the source terms describing a particle orbiting 
the star can be very easily accomplished. This is not the case for the Regge-Wheeler gauge, as even for the non-rotating case, 
one is forced to include second order derivatives of the source terms (Ruoff 2001b). Since the source terms contain 5- functions, 
one has to deal with second order derivatives thereof. In the axial case, no derivatives appear, and the perturbation equations 
with the source terms, which are given in Ruoff (2001b), are very simple. We expect the same to be the case for the polar 
equations in the BCL gauge. Here, it should be possible to plug the source terms into the equations for the extrinsic curvature 
without getting any derivatives. 

In subsequent papers, we will present results from the numerical evolution of the perturbation equations of slowly 
rotating relativistic stars in the BCL gauge with the particular focus on oscillations modes which are unstable with respect 
to gravitational radiation. As a step further, we will include the contribution of a test particle acting as a source of excitation 
for the stellar oscillations. 
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APPENDIX: THE PERTURBATION EQUATIONS FOLLOWING FROM EINSTEIN'S EQUATIONS 

Kojima (1992) derived the perturbation equations in the Regge- Wheeler gauge directly from the linearized Einstein equations 
without resorting to the ADM formalism. In this section we repeat this calculation using the BCL gauge. In order to facilitate 
the comparison with Kojima's equations, who uses the more familiar notation of Regge & Wheeler (1957), we switch to a 
similar notation. In the Regge- Wheeler gauge, the quantities ho and hi denote the axial perturbations of /it{0,e} and fer{(^>,e}) 
respectively, whereas the corresponding polar perturbations are set to zero. Since in the BCL gauge, the latter do not vanish, 
we denote them by /io,p and /ii,p, respectively, and to avoid confusion we denote the axial ones by ho,a and The remaining 
non-zero polar perturbations are then Hi and H2- Thus, the expansion of the metric in BCL gauge reads in this notation: 



/ -2u; {h'j^pd^, + h'oZsinede) 



Im 



V 



Hi" 

■k 
■k 



htr>do - hl^J sin 63^ h\]ld^ + h^";"^ si 



h'{%de - /i'Ca/ sin ed^ h'{%d^ + ft'i™ sin Ode 





(88) 



Note that here htt is not zero, which is a consequence of the relation between the perturbation of the lapse a and htt given 
by Eq. (^^. The relation between the above variables and the ones used in the previous sections is the following (we again 
omit the indices / and m): 



Hi = e^^Ki 

H2 = S3 , 

ho,p ~ Ks , 

ho,a ~ Kq , 

hi,p = V3 , 

hl,a = 14 , 

R = —ui , 

V = -U2, 

U = -U3. 

The extrinsic curvature components can be expressed as 

Ki = 2e-'''' {Hi-\'Hi+Lu[im{hlp^\'hi,p) +Ct^ {hl,-\'hi,a)]} 

r 2 

K2 — e~ h'o^p ho,p + Hi — hi^p + uj£^ /ii,a 



H2 — imujH2 



K-i = e ^"^ h'o^^ /io,a + hi,a - imujhi 



(89) 
(90) 
(91) 
(92) 
(93) 
(94) 
(95) 
(96) 
(97) 

(98) 
(99) 

(100) 



A very often occurring combination of variables in the perturbation equations is kg — hi for both the axial and polar cases, 
which we abbreviate with the following functions 



— ^0,a ~ hi, a j 



■0,p 



hi.p 



The equations coming from the {tt), (tr), (rr) and the addition of the (96) and (00) components can be written as 



42 + -^^ci 

with 



A-4 -^Im 



(I) 



±1 aW 



0, 



Im-^l — Im 



Z + lm Ai+i 



and 



^(") 



2e^ 
r2 

0, 

2 

+ - 



rH2 - Ah'i^p - IGnr^e^^C'^ {H ~ ^) + A (^X' - hi.p + ( 1 - 2t-A' + 



H2 



V 



C 



(tt) 



2uj (^\' 
-uj' +u(2u' + 2A' 

2A 



Za - 
2\A 



4lj 



h'o^a — 327rf2 (p + e) e 



21/+2A 



U 



(101) 
(102) 

(103) 

(104) 

(105) 
(106) 

(107) 



2co [Z'p -H[+ e'^H2) + {u;' - 2lu (A' + v')) h'o,p 



2loI u' + A' 



ho,; 
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uj' — 2uj { u' + \' 

r 



hi,p + {iu' + 2iu{X' - u')) Hi - 327106^"+^^ (p + e) V 



^H2 + ^{Zp + Hi- 2h'„,p + 2u'ho,p) + 167r (p + e) (e^'T? - Hi) 



2Au 



167rSl {p + e) 



hl,a , 



— + '^]H2~ 167vn {p + e) hi,p 



+ + \) '^^•^ - + ^) ^ - i: [(2-' + l) ^ f 



A 2A 

e 



rr 2l/ A / 

-"2 — e ^"•i.p I 



= 



ujh'o 



1^ H ^ I + 



-1/2 + 2e 



ii( + ( - - A' j ifi 
A 



2u-2\ / / 



ill 



A 



e hi 



A 



IdTve^^ip + e)H - (A + le^p) H2 + ^e'""'^ (i^' - A') 



= 2Lue~ 



^/io,a — Z'a + - A'^ (/lo.a — ^a)^ + 2uj'e ^"^ ^/lo,a " "^^o.a^ - IGTre^" t47(p + e)?7 , 
hip -Z'p + H[- e^"- (h2 + ^/io,p) + - ^') ('ii.p + ^1) 



+ e ^"^o;' (^Hi + — -/io.p^ — IG-kvjc^" {p + t)V 



The and {r0) components are 
Aa<i'+imd(i)+/:±Mi'+£±S<l)=0 



Im 

with 

(te) 



am 

;;(*«) 



~H2 + e- 



H[-Z'p + -h'o^p + (a' + u' ---'^Zp+ {u' - A') Hi- — (rX' 



rv + e 



1) ho,p 



+ 167re'" {p+e)V 



(108) 
(109) 
(110) 
(111) 

(112) 

(113) 
(114) 
(115) 

(116) 

(117) 
(118) 
(119) 

(120) 
(121) 

(122) 



2Aa; + - J^') hi,p) + ^' ( ^^^2 + 2/ii,p 



167rr2^x7(p + e) (ff2 + (l + C7')/f + C7''e) , (123) 



2uj'e ^^hi , 



V 

;;('■«) 



-Acje [h'i_a + {t/' - A') hi,a\ 



-Z: 

+ 



f)/^;„+[8.e-(p + e) + ^(l + rA')-^e^ 
//2 - 47rre^^ (p + e) -ff , 



ii,p 



167rr-^ti7 (p + e) {Hi + e^'i?) - ijA {Hi + Zp - ho,p) + cj' ( —H2 - 2 (A + 2) /lo.p ) , 
2uj'ho,a , 



From the {ttf)) and (r<j!>) components we get 

Abl2+im42+^'^:l+^\i2 = o 



with 



Z',+ (,.' + X'--)Za + -h', , + 16^6^"+^^ {p+,)U-\^(u' + X' 

\ r J r ir \ r 



.A-2 



ho,a 



g(t« ^ -SAcuh'i,^ + 



Alo hx' - u' - -) - {A - 2) LO' 



hi, a 



(124) 
(125) 



(126) 

(127) 

(128) 
(129) 

(130) 

(131) 
(132) 
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= 2loA (e'^J/a + {i^' - A') + lu' ( A/ii,^ - ^H!^ I + levrr^e^^ro (p + e) [H2 + (l + C,-^) H + C,"'^] 



2 2i^-2A,/ 21/ 

Za e hi a ~ e 

r 

Al^/lO,a + 2 (A+l)(^'- 

2aj'/io,p , 



A-2 2 _2A 

^:2~ + 7'^ 



r 



ho,a , 



^' ( A/io,p - ^^2 ) - Wnr^vo (p + e) {e^" R + Hi) 



From the (Off)) and the subtraction of (66) and (0<j!>) components we get 



Atim + imgim + C2 fir, 
with 







r I 2 -2A 

/ — uj r e 



Zp — —llD.p 



I 2 -2A 

g = —UJ r e 



Za ho.c 

r 



t 



— IGnr^zue^'^ (p + e)V 



> + £)U 



ho,p + {h'i,p + (ly' - A') 

2iy-2A 



-H2 — imujho,p 



ho, a + e^" (/i'l^a + (i^' - a') hi^a) — imu)ho,a 



(133) 
(134) 

(135) 

(136) 
(137) 
(138) 



(139) 
(140) 

(141) 
(142) 

(143) 
(144) 



These equations are fully equivalent to the ones derived within the ADM formalism. Although they still contain some second 
order derivatives, they can be easily brought by introduction of a few auxiliary variables into a first order hyperbolic form or 
even into characteristic form, which is very useful for the numerical evolution. 
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